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We study magnetic relaxation dynamics, memory and aging effects in interacting polydisperse 
antiferromagnetic NiO nanoparticles by solving a master equation using a two-state model. We 
investigate the effects of interactions using dipolar, Nearest-Neighbour Short-Range (NNSR) and 
Long-Range Mean-Field (LRMF) interactions. The magnetic relaxation of the nanoparticles in a 
time-dependent magnetic field has been studied using LRMF interaction. The size-dependent effects 
are suppressed in the ac-susceptibility, as the frequency is increased. We find that the memory dip, 
that quantifies the memory effect is about the same as that of non-interacting nanoparticles for the 
NNSR case. There is a stronger memory-dip for LRMF, and a weaker memory-dip for the dipolar 
interactions. We have also shown a memory effect in the Zero-field-cooled magnetization for the 
dipolar case, a signature of glassy behaviour, from Monte-Carlo studies. 
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I. INTRODUCTION 

Magnetism in nanoparticles has received enormous at- 
tention in recent years due to its technologicaliri^ as well 
as fundamental research aspects.—"—. Amid many stud- 
ies related to magnetic nanoparticles, the important con- 
cerns were related to the relaxation behaviour of the as- 
semblies of nanoparticles which has been addressed in 
the recent times ji"— The dynamics of an assembly of 
nanoparticles at low temperatures has gained a lot of at- 
tention over the last few years. In a dilute system of 
nanoparticles, the interparticle interaction is very small 
as compared to the anisotropy energy of the individual 
particles. These isolated particles follow the dynamics in 
accordance with the Neel-Brown modelj^ and the system 
is known as superparamagnetic. The giant spin moment 
of the nanoparticles thermally fluctuates between their 
easy directions at high temperatures. As the temperature 
is lowered towards a blocking temperature, the relaxation 
time becomes equal to the measuring time and the super- 
spin moments freeze along one of their easy directions. As 
the role of interparticle interaction becomes significant, 
the nanoparticles do not behave like individual particles; 
rather, their dynamics is governed by the collective be- 
haviour of the particles, like in a spin glassier— 

In our recent work,^^ we studied the slow dynamics of 
NiO nanoparticles distributed sparsely so that particle- 
particle interactions could be neglected. However, in re- 
ality, interparticle interactions play a major role in de- 
scribing various interesting phenomena observed experi- 
mentally in a collection of nanoparticles. This leads us 
to include interactions among the particles in our study. 
If the interparticle interaction is too small, the dynam- 
ics of an assembly of nanoparticles is a result of the in- 
dividual nanoparticle dynamics. Increasing the particle- 
particle interaction, the behaviour of the system becomes 
more complicated, even if we consider each nanoparticle 
as a single giant-spin. The dynamics of the assembly of 
nanoparticles involves various degrees of freedom related 
to the individual particles coupled to each other with in- 



terparticle interactions. In this paper, we will discuss 
various interactions comparatively. 

The main type of magnetic interactions that can play 
an important role in the assemblies of nanoparticles are: 
(i) long-range interactions (for example, dipolar interac- 
tions) among the particles, and (ii) short-range exchange- 
interactions arising due to the surface spins of the parti- 
cles which are in close contact. The dipolar interaction 
is a peculiar kind of interaction which may favour ferro- 
magnetic or antiferromagnetic alignment of the magnetic 
moments depending upon the geometry of the system. 
This interaction among the nanoparticles may give rise 
to a collective behaviour which may lead to a dynamics 
similar to that of a spin-glass. A resemblance from spin 
glass system owes to a random distribution of easy axes, 
which causes disorder and frustration of the magnetic in- 
teractions. The complex interplay between the disorder 
and frustration determines the state of the system and its 
dynamical properties. These systems are widely known 
as superspin glasses. This superspin glass phase has 
been characterised by observations of a critical slowing- 
dowufi^ a divergence in the non-linear susceptibility;^ 
and aging and relaxation effects in the low-frequency ac 
susceptibility The Monte Carlo simulations on the sys- 
tem of an assembly of nanoparticles show aging'^^ and 
magnetic relaxation behaviour"^^ like in a spin glass. For 
an interacting assembly of nanoparticles, the aging and 
memory effects are the two aspects that have been stud- 
ied extensively in recent years, however, mostly in the 
case of ferromagnetic particles i^Ti^iii^iiSr— 

In this paper, we investigate the effects of interactions 
on a collection of a few antiferromagnetic nanoparticles, 
which has received relatively lesser attention. Recently, 
we have shown that for NiO nanoparticles, a combined 
effect of the surface-roughness effect and finite-size effects 
in the core magnetization leads to size-dependent fluctu- 
ations in net magnetic-moment.^^ These size-dependent 
fluctuations in the magnetization lead to a dynamics 
which is qualitatively different from that of ferromag- 
netic nanoparticles. We will study the dynamics of the 
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FIG. 1: A lognormal distribution of nanoparticles with sizes 
ranging from 1.3oo-5.3ao, where ao = 4.17A. The width of the 
distribution is 0.6. The inset displays the magnetic moment 
vs the particle size for the same size range. The magnetic 
moment has a non-monotonic and oscillatory dependence on 
R. 



system by analytically solving the master equation for 
a two-state model. We will invoke various interactions 
and study the relaxation phenomena under these in- 
teractions comparatively. We will show the effect of 
the size-dependent magnetization fluctuations on the 
time-dependent properties of an interacting assembly of 
nanoparticles. We will also compare the dynamics with 
that of the non-interacting case. We will discuss the 
memory effects in the FC as well as the ZFC magne- 
tization protocols. The organisation of this paper is as 
follows. In Section |TT1 we discuss the two-state model. 
We discuss the ZFC and FC magnetizations in Section 
mil and the ac-susceptibility study in Section HVl In Sec- 
tion |Vl we show the memory effect investigations. The 
aging effect has been presented in Section I VII Finally, 
we summarise in Section [VIII 



II. RELAXATION PHENOMENA AND 
VARIOUS INTERACTIONS 

We assume a nanoparticle to be a giant spin in which 
the exchange interactions are so strong that all the spins 
in the particle show a coherent rotation in unison. Hence, 
the dynamics of a system of nanoparticles can be eas- 
ily understood as the dynamics of a giant spin-moment. 
Here, we are dealing with antiferromagnetic nanoparti- 
cles, which display size dependent fluctuations in the 
magnetization as shown in the inset of the figure FIG. 
[TJ On the other hand, the net magnetic moment of fer- 
romagnetic nanoparticles shows a linear dependence on 
the volume of the particlesj^ As we study the dynamics 
of antiferromagnetic nanoparticles, we might expect the 



role of these size-dependent fluctuations in the magneti- 
zation to be manifested in the time-dependent properties 
of these antiferromagnetic nanoparticles. We use a simple 
model where the energy of each particle i is contributed 
by the anisotropy energy and the Zeeman energy. Now, 
in order to incorporate interactions in the present model, 
we add an interaction energy term Vij in the Hamilto- 
nian. Thus, the Hamiltonian is written as: 
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and the energy of a particle is given by 
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where fli is the magnetic moment of the i*'^ particle, /x; = 
K, is the anisotropy constant, and H is the applied 
field. Because of the interaction term V^ , the energy 
of the i"^ particle depends on the state of all the other 
particles. We assume each nanoparticle to be an ising 
spin, and the magnetic moment of each particle to be 
fli = fJ-iSiZ, with Si — ±1. The various interactions used 
in the present study are: 

i) A Nearest-Neighbour Short-Range (NNSR) type 
interaction, which can be written as: 



(3) 



where the summation is taken over the nearest- 
neighbour sites. Here J^'' is the NNSR interaction 
strength. 

ii) A long-Range, Mean-Field (LRMF) type 
interaction,"^* which is given as: 
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where TV is the total number of nanoparticles and 
j'"' is interaction parameter. The above equation 
Eq. (jlj is scaled by N in order to insure that the 
total interaction-energy is proportional to N and 
not N^. 

iii) A long-range dipolar interaction, with the particle 
j separated by a distance from particle i, given 
as: 
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where a — ^ofi^/AirrQ characterises the strength of 
the dipolar interaction and rg is the lattice param- 
eter of the lattice, where nanoparticles are placed. 



3 



Variations in rp may result the nanoparticles to be 
getting closer or away from each other and hence, 
increasing or decreasing the effects of dipolar inter- 
actions. 

Let us consider a model case of NNSR. Using the mean- 
field approximation, we can find a mean field felt by the 
i^^ particle due to interactions with its neighbours. We 
can replace Eq. (jSj by its mean-field form as: 



E 

MF 
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where Hf^^ is the mean field on the i^^ spin from all 
the nearest neighbour spins. Thus, in this framework, 
each particle can be viewed under the influence of a local 
magnetic-field: 



Hf = H 
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The effect of interactions on a nanoparticle can be in- 
corporated by replacing the external field by an effective 
field which may lead to a mean-field equation of state 
and can be solved self-consistently. Hereafter, in all the 
cases of interactions, we will define an effective field 
on a nanoparticle i as the sum of the external field % and 
a locally-changing interaction-field i.e.. 
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The above equation shows that the role of the interac- 
tions is to modify the energy barrier, which is solely due 
to the anisotropy contributions of each particle in the 
non-interacting case. For strong interactions, their effects 
become dominant and the individual energy-barriers can 
no longer be considered to be the only relevant energy- 
scale. In this case, the relaxation is governed by a co- 
operative phenomenon of the system. The energy land- 
scape with a complex hierarchy of local minima is similar 
to that of a spin glass. We should note that in contrast 
with the static energy barrier distribution arising only 
from the anisotropy contribution, the reversal of one par- 
ticle moment may change the energy barriers of the as- 
sembly, even in the weak-interaction limit. Therefore, 
the energy-barrier distribution gets modified as the mag- 
netization relaxes. 

By defining an effective field given by Eq. (|8l), we can 
map an interacting assembly of nanoparticles as a collec- 
tion of non-interacting nanoparticles, where each particle 
experiencesan effective field H^^ which gets modified as 
the temperature changes. In the absence of the magnetic 
field, the superparamagnetic relaxation-time for the ther- 
mal activation over the energy barrier ICVi is given by 
Ti = To exp(/CT^i/fcBT), where tq, the microscopic time, 
is of the order of 10~^ s. The anisotropy constant K, has 
a typical value of about 4 x 10~^ J cm"-^ for NiO^^ 



The occupation probabilities with the magnetic mo- 
ments parallel and antiparallel to the magnetic field di- 
rection are denoted hy pi{t) andp2(i) = 1— respec- 
tively. The magnetic moment of each particle is supposed 
to occupy one of the two available states with energies 
-ICVi + filHf or -/CV- - n^Hf, where /if is the satu- 
ration magnetic-moment of the i"^ particle. These prob- 
abilities satisfy a master equation, which is given by. 



-pi(i) = Pi{t) + — 
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The magnetic moment of each particle of volume Vi can 
be written as fii = [2pi — 1] /if . as 



dt 
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We can define magnetization of each particle as Mi = 
fJ-i/Vi. Since Hf^ also contains a summation over all the 
other /ij's interacting with the particle, the right-hand 
side of the above equation can be written as the sum 
of the interaction term with j^^ particle and the term 
containing the magnetic field. Thus, we can write the 
Eq. (Uni) as: 



(11) 



Here A{i,j) represents the interaction term which is 
given below for each case separately, and aoi — 
~ (/if)' H/V,nT. 

(i) For NNSR interaction, A{i,j) is given as: 



-7-, ifi = j; 

Y-r T ^r,r+e7 ^ 7~ J 



where r is the position vector of any site and e is 
the nearest-neighbour distance. 

(ii) For LRMF interaction, we can write A{i,j) as 



2N ViTiT' ' J' 



(iii) For dipolar interaction we can write. 



There are N equations of the form given by Eq. pT|) 
for the system of N particles. We can write these N 
equations in a matrix form as: 



dt 



M = AM - ao, 
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where ^ is a A'' x iV matrix whose elements are given 
by A{i,j), defined above, and Ai and Og are vectors of 
length iV. The elements of Ai are the M^'s, and those of 
Oq are the agi's. Eq. ([T^ can be solved for the various 
interactions defined above. The formal solution of the 
matrix equation Eq. (jl2l) takes the form: 



M{t) 



-At 



C + O, 



(13) 



where {> = A ^ao, a time- independent solution of Eq. 
(ini) and C = 7W(0) - D, M{0) being the value of M at 
t — 0. Thus we can write, 



(14) 



The matrix A in the above equation Eq. (|14p is a general 
matrix of order N x N. To evaluate e-^*, we use the 
MATLAB function expm which implements a diagonal 
Fade approximation of exponential of the matrix with a 
scaling and squaring technique.—"^ 

Fade approximation to an exponential of a matrix B 
is written as: 



where the numerator term is given as 



{p + qV-jKp - jV- 



and the denominator term is written as 
[p + q- 3)W- 
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Diagonal Fade approximants can be calculated by taking 
q — pin Eqs. (fT6|) & (fT7| . The index p depends upon the 
norm of the matrix B. Using a backward-error analysis, 
p can be calculatedJ^ Unfortunately, the Fade approx- 
imants are accurate only near the origin. This problem 
can be overcome by the scaling and squaring technique 
which exploits, 



(18) 



Firstly, a sufficiently large k is chosen such that B/2'^ is 
close to the zero matrix. Then, a diagonal Fade approxi- 
mant is used to calculate exp(S/2'^). Finally, the result is 
squared k times to obtain the approximation to e^. We 
can use Fade approximation with a scaling and squar- 
ing technique to calculate e"^* in Eq. (fT4| . However, the 
norm of the matrix A, \\A\\ ~ 10^, is too large in all 
the three cases of interactions in Eqs. dSj, (Hj) & ([5]), for 
N = 1000. The scaling and squaring method comes very 
handy to overcome this large value of the normj^ In this 
case, the degree of the approximant is p = 13. The to- 
tal magnetic moment of the system of nanoparticles with 
volume distribution P{Vi) is given by 
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1 ^ r 



M,P{V,)dV,. 
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The size-distribution plays a significant role in the 
overall dynamics of the system of nanoparticles governed 
by equation (fT9|) . The exponential dependence of r on 
the particle size Vi assures that even a weak polydisper- 
sity may lead to a broad distribution of relaxation times, 
which gives rise to an interesting slow dynamics. For 
a dc measurement, if relaxation time coincides with the 
measurement time scale Tm, we can define^ a critical 
volume Vb as ICVb = A:B2Bln(Ti„/ro), where Tb is re- 
ferred as blocking temperature. The critical volume Ve 
has a strong linear dependence on Tb and weakly log- 
arithmic dependence on the observation time scale Tm- 
If the volume of the particle Vi in a polydisperse sys- 
tem is less than Ve, the super spin would have under- 
gone many rotations within the measurement time scale 
with an average magnetic moment zero. These particles 
are termed as superparamagnetic particles. On the other 
hand if Vi > Vb, the super spins can not completely 
rotate within the measurement time window and show 
blocked or frozen behavior. However, the particles having 
volume Vi ~ Vb are in dynamically active regime. The 
systems of magnetic nanoparticles are in general polydis- 
perse. The shape and size of the particles are not well- 
known but the particle-size distribution is often found 
to be lognormal4^ We consider the system consisting of 
lognormally-distributed and the volume Vi of each parti- 
cle is obtained from a lognormal distribution 
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where v — IniV), V is the mean size and a the width of 
the distribution. The distribution consists of 1000 par- 
ticles of sizes between R = 1.3ao and R = 5.3ao, where 
ao(= 4.I7A) is the lattice parameter of NiO*^ The total 
number of particles is deliberately chosen to be small in 
order to simplify the numerical calculation, as well as to 
get some insight into the role of size-dependent fluctua- 
tions in the magnetization in the present study. In all 
the cases, we use cr = 0.6 and v = 3.57. For the dipolar 
case, the nanoparticles are arranged inalOxlOxlO 
simple-cubic box. We can solve Eqs. p4)) and p9| for 
any heating/cooling process using any of the interactions 
discussed above. A recipe to solve these equations for a 
zero-field-cooled (ZFC) protocol is as follows. The sys- 
tem is cooled from a very high temperature to the lowest 
temperature in the absence of any magnetic field. The 
absence of the field during cooling causes a complete de- 
magnetization of the nanoparticles. This condition is the 
same as pi(0) = 1/2, or 7W(0) = in Eq. Now, a 

constant field is applied and the system is heated upto 
a high temperature. At each temperature change, we 
evolve the system using Eq. (fT4)) . Suppose we increase 
the temperature from T to T -|- dT . Then, the magneti- 
zation at T -f dT is given as: 



MzFc{t,T + dT) = e-^(^+''^)*7WzFc(0, T-fdT) 



(l-e-^^^+'^^^'MT + dT), 



(21) 
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FIG. 3: A comparison of the ZFC curves for all the interacting 
cases. Here, again a constant temperature decrease/increase 
of AT = 0.004 for every 100 s and a field of /i = 0.01. The 
values of interaction parameters are a — 2.5 x 10~®, J^'^ = 
2.5 X 10"^ and = 4.5 x 10"* for dipolar, NNSR and LRMF, 
respectively. 

A4 on T, which varies during the process. In the above 
expression, tj^ is the wait-time at each temperature- 
change. The magnetization is averaged over the volume 
distribution P{Vi) as: 

MzFC- J AhFc{t,T;V,)PiV,)dV,, (23) 

where MzFcit,T;Vi) is the i^^ element of MzFGit,T ~ 
dT), given by Eq. (Hi]). We can define the heating rate of 
the process as the total time elapsed at each temperature- 
change. Thus, a heating rate of 100 s per temperature- 
unit corresponds to a heating process in which the system 
is relaxed for t = 100 s at each temperature-step. 

Henceforth in this paper, volume Vi is used in units of 
the average volume V, which is IQSoq in our case. Also, 
the average anisotropic-energy ]CV is taken as the unit 
of energy. By setting /cb = Ij we can use /CF as a unit 
of temperature T and the field T-l. Hereafter, we use a 
dimensionless quantity h ~ ^iBH/K,a^ as the unit of the 
field, e.g., h — 0.01 is equivalent to a magnetic field of 
300 Gauss. 



FIG. 2: ZFC and FC curves for the (a) NNSR, (b) LRMF 
and (c) dipolar cases. For each case, there is a constant tem- 
perature decrease/increase of AT = 0.004 for every 100 s and 
a field h = 0.01. 

where the initial value of the magnetization vector at 
T + dT is the same as the magnetization vector at T, 
relaxed for the wait-time Thus, we can write: 

^ZFc(0,T-f dT) =7WzFc(i»,T). (22) 

Since the elements of the matrix A and the vector v de- 
pend on T, we have shown the functional dependence of 



III. ZFC AND FC MAGNETIZATIONS 

we begin our study with the most fundamental proto- 
cols, i.e., the study of zero-field-cooled (ZFC) and field- 
cooled (FC) magnetizations. In the ZFC process, the 
system is first demagnetized at a very high temperature 
and then cooled down to a low temperature in a zero 
magnetic-field. A small magnetic-field is then applied 
and the magnetization is calculated as a function of the 
temperature. During this heating process, the evolution 
of magnetization is given by Eqs (PTjl & (1^^ . On the 
other hand in the FC protocol, the system is cooled in 
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the presence of a small magnetic-field from higher tem- 
peratures to a low temperature. For a decrease of the 
temperature from T to T — dT , we can write: 

MFG{t,T-dT) = eA'^-'''^'>'MFc{0,T - dT) 
+ (l-e-^(^-'^^)*)j>(r-dT), 

(24) 

where 

MFciO,T~dT)^MFc{t^u,T). (25) 

The magnetization is averaged over the volume distribu- 
tion P{Vi), given as, 

Mfc = J MFc{t,T;V,)P{V,)dV,, (26) 

where Mpcit, T; Vi) is the i^^ element oi MFc{t, T~dT), 
given by Eq. ((M)) . In the present study, we calculate the 
ZFC and PC magnetization for various temperature, us- 
ing long-range and short-range interactions. For all these 
cases, we use a constant temperature decrease/increase of 
AT = 0.004 for every 100 s. In FIG.[l](a), (b) & (c), we 
have shown ZFC and FC magnetization against temper- 
ature plots for NNSR, LRMF and dipolar interactions. 
Since the interaction plays a role in decreasing the net 
magnetization,^'^ we have invoked an antiferromagnetic 
interaction in the case of LRMF and NNSR interactions. 
We find that on increasing the temperature, the magne- 
tization MzFC first increases, attains a maximum at a 
blocking temperature, and then starts decreasing. 

In all the cases, we find that as we increase the in- 
teraction strength, the peak of the ZFC magnetization 
shifts toward lower temperatures. In the case of LRMF, 
we find a dramatic smoothness in the ZFC magnetiza- 
tion as the interaction is increased. The smoothness of 
the ZFC magnetization indicates a lowering of the effects 
of size-dependent fluctuations in the magnetization. The 
reason for the smoothness may be the onset of a collec- 
tive behaviour of the nanoparticles with the increase in 
the interaction parameter. An increase in the interaction 
strength is equivalent to bringing nanoparticles closer to 
each other. As the nanoparticles come closer to each 
other, their dynamics loses their individuality and the 
net effect of other nanoparticles leads to a collective be- 
haviour. The same behaviour can be seen in all the other 
cases. For the sake of comparison, we have also shown 
the ZFC magnetization for the non-interacting case as 
well as for the interacting cases together in FIG. [H We 
find that in all the cases in FIG. [U the FC magnetiza- 
tion Mfc coincides with Mzfc higher temperatures, 
but departs from the ZFC curve at lower temperatures 
that are well above the blocking temperature. On further 
lowering the temperature, the FC magnetization tends 
to a constant value. The blocking temperature shows 
a substantial dependence on the heating rate. For an 
infinitely-slow heating-rate, Tb approaches zero and the 
ZFC curve shows a similar behaviour as the FC curve. 



We also find that for a very weak interaction, the FC 
magnetization never decreases as the temperature is low- 
ered. But, the increase in the interaction leads to a flat- 
ness in the FC magnetization below the blocking temper- 
ature which can be seen in all the cases shown in FIG. 
[2] This flatness in the FC-curves again shows the signa- 
ture of the co-operative phenomenon of the nanoparti- 
cles, arising due to the interactions. We see that the FC 
and ZFC magnetizations in the non-interacting case are 
greater than those in the interacting cases. We also find 
that the blocking temperature for the non-interacting 
case {Tb = 0.04) is much lower than the corresponding 
interacting one {Tb = 0.05). 

IV. AC-SUSCEPTIBILITY 

The relaxation phenomena of a magnetic system are 
often investigated in the presence of an oscillating 
magnetic-fieldi^ii^'^i^'^ It is worthwhile to use an oscil- 
latory magnetic-field in our study as the frequency of os- 
cillations of the magnetic field may compare with the in- 
verse of the relaxation time of the nanoparticles. Hence, 
by changing the frequency of the magnetic field, we can 
study the dynamics of an assembly of nanoparticles. For 
a time varying magnetic field h{t) — /loe"*"^*, the ana- 
logue of Eq. is given by: 

^M=AM-d{t), (27) 
dt 

where d{t) = aoe~*'^*, and w is the frequency of oscillation 
of the magnetic field. Multiplying by e"-^* on both the 
sides of the above equation Eq. ([?7)) gives: 

_ AM^ - -e-^*d{t), (28) 

which can be written in a simpler form as: 

^ (e--^*^) = -e--^'d{t) (29) 

The solution of the above differential equation is given 
by, 

M = e-^'C ~ f e^(*-*)a(s)ds, (30) 
Jo 

where C — TWq. Using d{s) = aoe~*"^ in the above 
equation, we get, 

M = e^'Ma - I e^^'-'^aoe-'^'ds. (31) 
Jo 

Carrying out the integral in the above equation, we have, 

M = e-^'Mo + {A + iujir^doe-''^' 

- e^\A + iujI)-^do, (32) 
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FIG. 4: A comparison of tlie (a) real, and (b) imaginary com- 
ponents of tlie ac-susceptibility for tiie non-interacting and 
interacting cases for botii the cases, tlie frequency of the ac- 
field is cj/27r = 1/500. 
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FIG. 5: (a) Real and (b) imaginary parts of the susceptibility 
for the interacting case (LRMF) are plotted against the tem- 
perature for various frequencies of the applied ac-magnetic- 
field. 



where I is the identity niatrix of order N x N. Eq. ([5^ 
involves the matrix exponential term e"^*, which can be 
evaluated using diagonal Pade approximants with a scal- 
ing and squaring technique from MATLAB. This proce- 
dure has been discussed in Section [TTl Here, using matrix 
exponential function expm, we can calculate e"^* in Eq. 
([5^ . From Eq. (HH), we can easily calculate M{t) by 
averaging the above value of the magnetization over a 
size-distribution. The response of the time- varying field 
h{t) = /ioe~*"* can be expressed as the sum of an in-phase 
component, x' ^^^d an out-of-phase component, x", of 
the susceptibility x(a;) = x' ~ *x"- The in-phase compo- 
nent with the applied field, x', is the lossless component, 
while the out-of-phase component with the applied field, 
x", is the lossy component. It is obvious that if the 
change of the external field is very fast as compared with 
the relaxation time of the particles t, (w ^ t~^), then 
the particles cannot follow the field variation; hence, the 
magnetization gradually decreases to zero as uj increases. 
This term is maximum at a; = and decreases gradu- 



ally as w increases or decreases from the point uj ~ t^^. 
If the frequency is much less than than re-orientation of 
the magnetization (i. e. w ^ t~^, then the magneti- 
zation is always in equilibrium over the time-scale of the 
measurement. The dynamic susceptibility can be defined 
as: 



Xit) = M{t)/h{t). (33) 
The Fourier transform of x(i) can be written as: 

XH ^ x'(^) - *X"(^) ^jj^ X{t)dt. (34) 

Now using M{t) = Mn + iMj, where and Mj are 
the real and imaginary parts of the magnetization, we 
can separate x' and x"> respectively, as: 



1 



T 



X'(w)-r-^/ (AfflSin(wt) -KMicos(a;t)), (35) 

' Jo 
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and 

1 r 

= — — / (AfflCOs(a;t) - Misin(a;t)). (36) 
hoT Jo 

Using the above set of equations, we can calculate the 
real and imaginary parts of the ac-susceptibility. For 
the ac-susceptibility calculations, the system is demagne- 
tized at the lowest temperature and an ac-magnetic-field 
of amplitude ho = 0.01 and frequency uj/2tt = 0.01 is 
applied. The system is then heated upto a high temper- 
ature. At each temperature change, we evolve the system 
for T = 100 s and calculate x' and x" using the above dy- 
namical equations. It has been experimentally observed 
that ac-susceptibility measurements are sensitive to the 
interaction effectsi^ and confirmed using Monte Carlo 
simulations.'^^ We also compare the effect of interactions 
by considering the non-interacting and interacting cases 
individually. We find that non-interacting case, which is 
well described by the individual relaxation of particles 
differs qualitatively from the interacting case, where the 
interactions among the particles make the dynamics com- 
plex. The real and imaginary parts of ac-susceptibility 
for the non-interacting case as well as the interacting 
case using LRMF interactions is shown in FIG. HI For 
both the cases, the frequency of the ac-field is fixed to 
u}/2tt = 1/500. We find that the due to the interaction, 
the in-phase part as well as the out-of-phase part show 
reduction in the net value of the susceptibility. The ef- 
fect of size-dependent fluctuations can be observed in the 
in-phase and out-of-phase components of the susceptibil- 
ity. The effect of these size-dependent fluctuations are 
displayed in the non-interacting as well as the interact- 
ing cases, as ripples in the ac-susceptibility. But, it is 
interesting to see that the ripples are smoothened in the 
interacting cases. The correlated behaviour of nanopar- 
ticles is responsible for the averaging of the ripples in this 
case. The effects of interaction can also be seen as the 
shifting of the peaks of x' and x" curves towards lower 
temperatures. 

We have also performed a detailed study of the fre- 
quency dependence of the real and imaginary parts of 
the susceptibility. In FIGs. \5\ (a) and (b) we have shown 
a frequency dependence of the susceptibility for the in- 
teracting case by plotting x' and x" versus temperature 
at various frequencies. We see that the increase in the 
frequency causes a lowering in the magnitude of x' and 
x" ■ Again, we find that the ripples in the susceptibility 
components due to size-dependent fluctuations get sup- 
pressed with an increase in the frequency. We also see 
that the peak in the x' and x" curves slightly shifts to- 
wards higher temperatures as the frequency is increased. 
This frequency dependence shows that the particles be- 
come less responsive to the magnetic field as the fre- 
quency is increased. For a very large frequency the par- 
ticles may not flip due to the magnetic field which may 
result in zero-magnetization. 



V. MEMORY EFFECTS 

A. FC memory-efFect 

Sun et ofi^ have shown a memory-effect phenomenon 
in the dc-magnetization by a series of measurements 
on a permalloy NigiFeig nanoparticle sample. Later, 
this memory effect was also been reported by Sasaki 
et alM> and Tsoi et al^ for the non-interacting or 
weakly-interacting superparamagnetic system of ferritin 
nanoparticles and Fe203 nanoparticles, respectively. All 
of the studies were concerned only to the ferromagnetic 
nanoparticles. Experiments on NiO nanoparticles by 
Bisht and Rajeev^^ also confirm a weak memory effect 
in these particles. Chakraverty et al^ have investigated 
the effect of polydispersity and interactions among the 
particles in an assembly of nickel ferrite nanoparticles em- 
bedded in a host non-magnetic Si02 matrix. They found 
that either tuning the interparticle interaction or tailor- 
ing the particle-size distribution in a nano-sized magnetic 
system leads to important applications in memory de- 
vices. In our recent study,'^'* we had also shown the 
memory effects for an assembly of non-interacting NiO 
antiferromagnetic nanoparticles by the analytical solu- 
tion of a two-state model. The protocol for the memory 
effect is as follows. Firstly, we cool the system from a 
very high temperature to Tbaso = 0.005 with a probing 
field of h = 0.01 switched on. The system is again heated 
from Tbaso to get the reference curves, which are shown 
as "Ref." in FIG. [H The dynamics of the system in this 
process is the same as that given by Eqs. (|^ & 

during the coohng process, and Eqs. (PT|) . (1^ & 
during the heating process. During these processes, 
we allow 100 s to ellapse for every temperature step of 
AT — 0.001. We again cool the system from Th to Tbaso 
but this time with a stop of 10000 s at T = 0.04. The field 
is cut off during the stop. At the stop, the magnetization 
is given by, 

^oooi(t, T,t,op) = e-^(^--)*>lcooi(0, Tstop), (37) 

where 

M cool (0,Tstop) = A^cool(i«;, Tstop). (38) 

and 

Afoooi = J M,ooiit,T;V^P{V,)dV,. (39) 

After the pause, the field is again applied and the system 
is again cooled upto the base temperature Tbaso- The 
recovered magnetization can be given as: 

^oooi(t, T,t,op - dT) = e-^(^---'*^)*7Wcooi(0, T,top - dT) 

+ (1 - e-^(^---^^)*)z>(T,top - dT), 

(40) 

where 

^oooi(0,T,top - dT) = Mcooi{tu„Tstop)- (41) 
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FIG. 6: PC memory-effects for various interacting cases and 
the non-interacting case with a stop of 10000 s at T = 0.04. 
The uppermost set of curves shows the non-interacting case, 
while the dipolar and the NNSR cases follow them towards 
the lower magnetization. The LRMF case is shown at the 
bottom. 



The process is shown as the "cooHng" curves in FIG. IH) 
Finally, we heat the system at the same rate as that of 
cooling without any stop, which is shown as the "heating" 
curves in FIG. [51 The heating process follows Eqs. (HU, 
((22)) & ((23l) . However, the transition from the cooling 
process to the heating process at Tbasc is given by. 



A^hcat(i,Tbaso + dT) 



:-^(^---^^)*)^(rba.c 



where 



M heat (0, Tbasc + dT) = heat Tbasc)- 



+ dT) 
dT), 
(42) 



(43) 



As has been discussed earliar that the interactions 
among the particles play an important role in the slow 
dynamics, we study the memory effect using the above 
protocol for various long-range and short-range interac- 
tions. Since the memory-effect phenomenon in the FC 
inagnetization protocol arises due to polydispersity in the 
system of nanoparticles, we find a memory dip in all the 
cases. Yet, we can find out how strong the memory effect 
is under various interactions. We compare the relative 
strength of the memory effects for various interactions by 
introducing a parameter, the Memory Fraction (M.F.), at 
the stop during the memory measurement, defined as: 



Memory Fraction (M.F.) 



AM M - AT 



NT 



NT 



(44) 



The calculated values of the M.F. for the no interaction, 
dipolar, NNSR and LRMF cases are given in Table ID In 
all the cases, we fix the waiting-time at the stop to be 
10000 s. The trend is a bit surprising. As we know that 
the memory effect in non-interacting antiferromagnetic 



Interaction 


M.F. 


Non-interacting 


0.24 


Dipolar 


0.18 


NNSR 


0.22 


LRMF 


0.302 



441). has 



TABLE I: The memory fraction, defined by Eq. 
been given for the various interactions. 



nanoparticles is more than that for ferromagnetic parti- 
cles of the same size-distribution when the distribution 
consists of smaller sizes.— Here, we find that invoking 
various interactions do change the net memory-dip in an- 
tiferromagnetic nanoparticles. As compared to the non- 
interacting case, we see a reduction in the memory dip for 
dipolar interactions, an almost similar dip for the NNSR 
case, and an increased memory dip for the LRMF case. 
Since we are employing a simple model in the present 
study, we find the best effects of the interactions in the 
case of LRMF. We find that as the interaction param- 
eter increases, the effect of the size-dependent fluctua- 
tions decreases, and thus, the memory dip decreases in 
the dipolar interactions case. On further increasing the 
interaction parameter, a collective dynamics becomes re- 
sponsible for the enhanced memory-dip in the NNSR and 
LRMF cases. 

In our earlier work,'^'* we had analysed the role of poly- 
dispersity in the memory effect for the non-interacting 
case. As there was no interaction, the dynamics of the 
individual particles was responsible for this peculiar ef- 
fect. However, in the case of interacting antiferromag- 
netic nanoparticles, we can see that the size-dependent 
magnetization-fluctuations and interactions among the 
nanoparticles do modify the memory dip. This memory 
effect in an interacting system of nanoparticles can be un- 
derstood in terms of the droplet picture.— Droplet theory 
earlier, proposed for spin-glass systems, is very helpful 
in getting some insight into spin configurations in these 
interacting systems of nanoparticles. According to this 
model for a spin-glass, if the system is rapidly quenched 
in a field /i to a temperature T below the critical tem- 
perature, spin-glass domains, or clusters called droplets, 
which are in local equilibrium with respect to T, and h 
grow in size. In this picture, a small temperature-change 
causes substantial changes in the equilibrium state. At 
any time t, droplets of various sizes exist. Thus, the 
system is analogous to the non-interacting case where 
various particle-sizes can be replaced by various clus- 
ters grown at temperature T and field h. In droplet 
theory^i^ the dynamics of the droplets is considered to 
be a thermally-activated process where the energy barrier 
to form a droplet scales with L, the size of the droplet. 
This droplet picture is relevant to the present study. 
The thermally-activated dynamics of droplets shows a 
resemblance to the two-state model of the superparam- 
agnet. We see that due to the mean- field interaction with 
all the other particles, the individual magnetic moment 
can no longer remains dissociated from the other spins. 
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FIG. 7: ZFC memory-effect using Monte Carlo simulations. 
The difference between the aged and the normal (reference) 
ZFC magnetization is shown. The memory dip around the 
stop temperature Ts = 0.1 can be seen. In the inset, we 
have shown the aged and the reference ZFC-magnetizations 
together. 

Hence, an effective moment replaces the individual mo- 
ment. Since the averaging is done due to the presence 
of other spins, the net effect on the dynamics is modified 
according to the interactions between them. As the mem- 
ory dip AM depends strongly on the wait-time param- 
eter and the distribution of sizes the parameter AM/Af 
becomes very effective in the qualitative description of 
the FC memory-effect. 



B. ZFC memory effect 

Sasaki et alM suggested a ZFC memory-effect protocol 
to confirm whether the observed memory effect is due to 
glassy behaviour or not. In this method, we first cool the 
system rapidly in a zero-field from a high temperature 
to the stop temperature, which is well below the block- 
ing/freezing temperature. The system is left to relax for 
a wait-time t^. The rapid cooling is then resumed down 
to the lowest temperature where a small magnetic field is 
applied, and the magnetization is calculated during the 
heating process. A reference ZFC-curve can also be ob- 
tained without any stop during the cooling process. If 
system exhibits memory effect, a memory dip should be 
observed around the stop temperature during the heat- 
ing process. Since this dip is substantially smaller, it is 
better to see the behaviour of the difference between the 
aged and the normal ZFC magnetizations as a function 
of the temperature. In the present context of study, we 
have also investigated the ZFC memory-effect. The sys- 
tem, during sudden cooling under a very small field of 
0.000001 (almost zero-field), is put on hold for a stop 
of 10'* s at T = 0.04. During heating the system, we 
find a very small dip at the stop temperature. This very 
small dip is a drawback of our simplified model. Since we 



did not incorporate randomness in our model, the simple 
two-state model could not capture the ZFC memory-dip. 
In order to demonstrate the memory effect, we perform 
Monte Carlo simulations^— in a 5 x 5 x 5 simulation box 
of a simple cubic lattice with 125 particles of different 
sizes. Again, the volume Vi of each particle is obtained 
from a lognormal distribution as defined in Eq. (|20|) . In 
the Monte Carlo method, during each Monte Carlo step 
(MCS), we select a particle i from 1000 random sequences 
of 125 particles with either 'up' or 'down' superspin ori- 
entations. The attempted fiip of the orientation is ac- 
cepted with a probability of exp(— iJi/fc^T), where Ei is 
the energy of the i*^ particle given by Eq. ([2]). A Monte 
Carlo steps is the time-unit in our simulation. Whenever 
a flip attempt is successful, the magnetic moment of the 
particle i is updated according to the magnetic-moment 
versus particle-size plot shown in the inset of FIG. [T] One 
Monte Carlo step is equivalent to the intrinsic time-step 
To of the superspins. All the processes discussed below 
using the Monte Carlo method are under fast cooling and 
heating rates. Though very fast cooling and heating rates 
are seemingly far from reality, a qualitative picture of the 
interacting nanoparticles evolves by this method. 

In the present study, we focus on the dipolar- 
interaction term and take the interaction parameter to 
be a — 5x 10^^. The variation of the interaction param- 
eter a has the same effect as that by changing the lattice 
constant of the lattice on which the particles are placed. 
We repeat the same protocol for the ZFC memory-effect 
as discussed above. We first cool the system rapidly in 
a zero-field from a high temperature T — 0.6 to the stop 
temperature Ts = 0.1 (well below the blocking/freezing 
temperature) and let the system relax for 50000 MCS. 
The rapid cooling is then resumed down to the lowest 
temperature where a magnetic field of h — 0.01 is ap- 
plied and the magnetization is calculated during the heat- 
ing process where the temperature is changed in steps of 
AT = O.OlKV/kB for every 100 MCS . A reference ZFC- 
curve is also obtained without any stop during the cool- 
ing process. We have plotted the difference between the 
aged and the normal ZFC magnetizations as a function 
of the temperature for the interacting case with dipolar 
interactions among the particles in FIG. [71 For the non- 
interacting case we pointed out no ZFC memory-effect 
for a non-interacting assembly of NiO nanoparticles^ 
however, for the present case we see a dip appearing 
in AM — M^ef — M just at the temperature T, = 0.1, 
where the stop and wait has been made during the ZFC 
rapid cooling process. The ZFC memory-effect can be ex- 
plained by droplet theory.^'^ It has been reported earlier 
that below a critical temperature droplets (clusters) are 
growing as the time elapses. This time-dependent growth 
of the clusters responsible for the unusual memory-effect 
in the interacting assembly of nanoparticles. As we allow 
more time to elapse at Tg, the clusters grow proportion- 
ally to the wait-time. When the cooling is resumed, this 
equilibrated clusters freeze at lower temperatures. Dur- 
ing the heating process, just near Ts, the frozen clusters 
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FIG. 8: A comparison of the aging effect for four cases: No 
interaction, Dipolar, NNSR, and LRMF. In each case, the 
wait-time is — 100 s. 



rearrange themselves, and a dip can be observed in the 
AM vs. T curve. 



VI. AGING EFFECTS 

The magnetization relaxation of nanoparticles can also 
be studied by analysing the wait-time dependence of 
magnetization. In one such protocol, the system is cooled 
from a high temperature to the lowest temperature in 
the presence of a constant magnetic-field. During the 
cooling process, the system is halted at a temperature 
lower than the blocking temperature. The system is 
paused for a wait-time tw before switching off the mag- 
netic field. The relaxation of magnetization show the 
dependence. This phenomenon is known as aging effect. 
Aging effects in polydispersed assemblies of interacting 
nanoparticles have been been receiving a great deal of 
attention. The aging effects in interacting nanoparticles 
has been indicated a glassy behaviour in these systems. 
The aging-effect has been well-studied in the spin-glass 



systems 



54-56 



as well as in interacting assemblies of fer- 



romagnetic nanoparticles The aging effect or 

wait-time dependence in the thermoremanant magneti- 
zation (TRM) protocol has also been observed in non- 
interacting assemblies of nanoparticles. Most of the aging 
effect studies are centred around ferromagnetic nanopar- 
ticles. However, recent experiments on antiferromagnetic 
NiO nanoparticles show aging effect in these systems 
also.''^ In our earlier work,'^'* we had shown the aging ef- 
fect for a collection of a few non- interacting NiO nanopar- 
ticles. We find that in the non-interacting case, the wait- 
time dependence was due to an ill-defined initial state of 
the system just before switching off the field. When the 
interaction among the particles comes into the picture, 
we should be able to observe the effects arising due to 
the collective behaviour of the nanoparticles. Most of the 
aging experiments are performed by the measurement of 



the time-dependent ZFC and TRM magnetizations. The 
methodology in a TRM protocol is as follows. The sys- 
tem is cooled in a field to a base temperature Tbase below 
the blocking temperature Tb. During the cooling process, 
the magnetization is given by Eqs. (HH), ([25]) & (dH). Af- 
ter a waiting-time ^w, the magnetic field is switched off 
and we find the relaxation in the magnetization. The 
magnetization in this case is given by. 



M 



TRM 



{t, Tstop) = e-^(^--)*>lTRM(0, Tstop), (45) 



where 



and 



A^TRm(0, Tstop) — A^TRM(itu, Tstop) 



Mtrm= / MTRMit,T;Vi)PiVi)dV,. 



(46) 



(47) 



when the TRM magnetization is plotted against the time, 
it is found that the relaxation in the magnetization de- 
pends upon the waiting-time t^. 

In this section we will study the aging effect for a 
polydispersed system of interacting NiO nanoparticles in 
the TRM protocol. Our investigation is carried out by 
cooling the system in the presence of magnetic field of 
h ~ 0.01 upto the base temperature Tbaso = 0.024 and 
then cutting the field off after a wait-time tw to let the 
system relax. In a previous work,'^'* we had discussed the 
aging effects for various size-distributions and made a 
comparative study with the ferromagnetic particles case. 
However, in this section, we perform a comparative study 
of the various interactions in NiO nanoparticles, as dis- 
cussed in Section ini In FIG. [51 we plot the thermorem- 
anant magnetization against the time on a logarithmic- 
scale for the non-interacting as well as the interacting 
cases: dipolar, NNSR and LRMF. Before switching off 
the field, we wait for — 100 s. We see that the de- 
cay in the magnetization is fastest in the non-interacting 
case, which shows a two-step relaxation.^'* In the non- 
interacting case, the dynamics of an individual particle 
is not correlated with that of others in its surroundings. 
The role of size-dependent fiuctuations can be seen as 
ripples in the relaxation of the magnetization, but in the 
interacting cases, we see a smoother curve which shows 
the averaging of the fluctuations in the magnetization 
due to the mean field arising from the neighbouring par- 
ticles. Due to the interactions the magnetization per- 
sists for a longer time as compared to the noninteracting 
case. Hence, we can say that due to the interactions, 
the dynamics of the nanoparticles gets slower. We have 
also plotted the magnetization versus logarithmic-time 
for various wait-times — 100 s, 1000 s, and 10000 s in 
FIG. O We see a wait-time dependence in all the cases, 
though weak. In the inset of these curves, we have shown 
the corresponding relaxation-rate for each case, defined 
as —dM/d\ogiQ{t). Multiple peaks in all the cases refiect 
the distinction between the antiferromagnetic nanopar- 
ticles and the ferromagnetic nanoparticles: an effect of 
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FIG. 9: The aging effect for waiting-times t-w = 100 s, 1000 s and 10000 s has been shown for: (a) no interaction, (b) dipolar c) 
NNSR, and d) LRMF interactions. In the interacting cases (b) (c) and (d), we see a smooth decay of the magnetization with 
time. 



size-dependent magnetization- fluctuations. In a nutshell, 
we conclude that the role of size-dependent fluctuations 
in the magnetization dynamics has been controlled a little 
bit due to the interactions. The decay of the thermore- 
manant magnetization gets slower as the interaction is 
increased. 



VII. CONCLUSIONS 

We have studied the slow dynamics of an interacting 
assembly of antiferromangetic nanoparticles by solving a 
two-state model analytically. A collection of a few in- 
teracting antiferromagnetic nanoparticles has been com- 
paratively studied using various long-range and short- 
range interactions. We find that the ZFC magnetization 
shows ripples in all the cases, which starts decreasing as 
the interaction parameter increases. Due to the interac- 
tions, the dynamics of the system is governed not only 
by a broad distribution of particle relaxation-times aris- 
ing from the polydispersity, but also by a collective be- 
haviour amongst the particles themselves. We have also 



studied the dynamics of the system using an ac magnetic- 
field. We find that the real and imaginary components of 
the ac-susceptibility show ripples, which are due to the 
size-dependent magnetization-fiuctuations. The effects 
of these size-dependent fiuctuations are displayed in the 
non-interacting as well as interacting cases; however, we 
find that the ripples are smoothened in the interacting 
cases. Due to the interaction, we also find a shifting of 
the peaks of the x' x" curves towards lower tem- 
peratures. We have also calculated the frequency de- 
pendence of x' s-nd x", and found that the frequency 
causes a lowering in the magnitudes of x' x" ■ The 
frequency dependence of x' and x"shows that the parti- 
cles become less responsive to the magnetic field as the 
frequency is increased, we have also shown the memory- 
effect in the field cooled magnetization protocol for an 
interacting polydispersed assembly of nanoparticles. We 
have studied the effects of interactions on the memory 
dip by varying the interaction parameter. As compared 
with non-interacting case, we find a reduced memory- 
dip for dipolar interactions, an almost same memory- 
dip for NNSR interactions, and an increased memory- 
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dip for LRMF interactfons. Due to the limitations in our 
model, we find the best effect of the interactions only in 
the case of LRMF. We find that as the interaction pa- 
rameter increases, the effect of the size-dependent fluc- 
tuations decreases, and thus, the memory dip decreases 
in the dipolar interactions cases. On further increasing 
the interaction parameter, a collective dynamics becomes 
responsible for the enhanced memory dip in the NNSR 
case and LRMF cases. The memory effect and a weak 
aging effect have also been discussed for all the cases. 
For the non-interacting case, the size-dependent fluctua- 
tions play a role in the obtaining of a high memory frac- 
tion. As the interactions are incorporated, the memory 
fraction decreases, but for the strong interaction cases. 



the memory fraction increases. As the present model 
could not incorporate randomness and disorder, we have 
performed a Monte Carlo study considering the dipo- 
lar interactions among the particles. We find a memory 
dip in the ZFC magnetization, which indicates the im- 
portance of the interactions among the particles. The 
ZFC memory-effect is attributed only to the collective 
dynamics of the nanoparticles. The system of interact- 
ing nanoparticles under various interactions also shows 
a wait-time dependence. The decay of thermoremanant 
magnetization slows down as the interaction among the 
nanoparticles increases. We have done a comparative 
study of the aging effects in the no interaction, dipolar, 
NNSR and LRMF cases. 
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